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ABSTRACT 


Improved parameters are presented for predicting 
the tropospheric effect on electromagnetic range measure- 
ments from surface meteorological data. More geographic 
locations have been added to the earlier list. Parameters 
are given for computing the dry component of the zenith 
radio range effect from surface pressure alone with an 
rms error of 1 to 2 mm, or the total range effect from the 
dry and wet components of the surface refractivity, N, 
and a two-part quartic profile model. (N = 10 6 (n - 1), 
where n is the index of refraction. ) The new parameters 
are obtained, as before, from meteorological balloon data 
but with improved procedures, including the conversion 
of the" geopotential heights" of the balloon data to actual 
or "geometric" heights before using the data. The re- 
vised values of the parameter k (dry component of verti- 
cal radio range effect per unit pressure at the surface) 
show more latitude variation than is accounted for by the 
variation of g, the acceleration of gravity. This excess 
variation of k indicates a small latitude variation in the 
mean "molecular weight" of air and yields information 
about the latitude -varying water vapor content of air. Sta- 
tistical values of the height integrals of both dry and wet 
components of N are given for all the locations studied. 
Some comparisons of the values of these integrals are 
shown for stations a few hundred miles apart. 
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1. INTRODUCTION AND BACKGROUND 


The lower, nonionized atmosphere (troposphere 
plus stratosphere) retards the passage of electromagnetic 
signals and thus introduces an error into any electromag- 
netic measurement of range or range rate. The observed 
range is the line integral, along the signal path, of the in- 
dex of refraction, n, which at sea level differs from unity 
by 3 or 4 parts in 10^ and which varies along the path. In 
an atmosphere that is radially symmetrical about the earth 
(no horizontal gradients of n), a signal propagated verti- 
cally undergoes no path bending, and the observed range 
is Jn dh, where h is the height. The range error caused 
by refraction is then J(n - 1) dh. 

Range errors for other directions of propagation are 
related to the error at the zenith, but at low angles they de- 
pend on the profile shape also. At any angle, however, 
knowledge of the effect at the zenith is one of the require- 
ments for predicting the correction for measurements ex- 
tending through the atmosphere (e.g. , to earth satellites 
and more distant objects). Much of the work reported here 
deals with predicting the correction at the zenith on the 
basis of local conditions at the earth's surface at the time 
of the range measurement. This work is a continuation of 
a long-term study. 

The effect of the nonionized atmosphere is produced 
mostly by the troposphere proper. For brevity, the term 
"tropospheric range error" will be used here to include the 
stratosphere also, unless otherwise indicated. 

The refractivity, N, of air (defined as l()6(n - 1)) 
can be conveniently expressed as the sum of two parts, the 
so-called "dry" and "wet" components, here subscripted 
"d" and "w. " In the radio range (up to at least 15 GHz), 
both components are significant and are independent of 
frequency. For optical signals, on the other hand, es- 
sentially the whole effect is due to the dry part, and there 
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is some wavelength dependence. The values of computed 
from the Smith and Weintraub radio expression (Ref. 1) 
and from the Edle'n optical expression (Ref. 2) are in very- 
good agreement in the near -infrared region (wavelengths of 
a few micrometers). Values in the visible region differ 
from these by a few percent. The Smith and Weintraub ex- 
pressions for the radio region will be used here. The work 
on the dry component, as reported below, can be applied 
to both wavelength regions, with a wavelength correction 
of a few percent for optical uses, or the addition of the wet 
component effect for radio uses. 

The Smith and Weintraub expressions (Ref. 1) for 
the radio refractivity of air may be written: 



( 1 ) 


where T is the Kelvin temperature, P is the total pres- 
sure, and e is the partial pressure of water vapor, the 
latter two in millibars. 

The refractivity components, N d and N w , decrease 
with height above the surface of the earth, but at different 
rates. Theoretically would decrease exponentially with 
height in an isothermal atmosphere. Exponential expres- 
sions frequently have been used to approximate the N pro- 
file, both the total N and the separate components (Refs. 

3 and 4). 

It was previously shown (Ref. 5) that, in a dry at- 
mosphere, if the air temperature varies linearly with 
height (constant lapse rate a, where a= - dT/dh), the 
N profile is theoretically a polynomial function of height, 
not exponential. The smaller the lapse rate, the higher 
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the degree of the polynomial (approaching an exponential 
as OL -* 0). If heights are measured above the surface, the 
general form of the equation is: 


N = 



h * h, 
d 


where 



( 2 ) 


(3) 


M 


= 

Ra 


1 , 


(4) 


and the subscript s refers to surface values. In Eq. (4) 
g is the acceleration of gravity and R is the gas con- 
stant per gram of dry air. 


If a = 6. 8°C/km, which is near its usual value 
in the troposphere, then H = 4 for the profile. Using 
the same form for the N w profile, for convenience, the 
two-part expression becomes (Ref. 5): 



N 

w 


N a 

rf 1 < h d - h > ■ 

h , 


N 

w 



w 



h s h , 
d 


h £ h 
w 


(5) 


It was also shown previously (Ref. 6), both theo- 
retically and from observed data, that although at any 
point is a function of both pressure and temperature, its 
zenith integral is a linear function of surface pressure 
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only. The tropospheric contribution to a vertical range 
measurement (electromagnetic) through the atmosphere is: 

Ah tro = 10' 6 J'N d dh = kP s , (6) 

where k is a constant for a given location. 

To reconcile Eqs. (5) and (6), it is necessary 
(Ref. 6) that 


where h d is height above the surface, a d is a constant, and 
the surface temperature, T , is in degrees Kelvin. This 
linear variation is also implied in Eq. (3); thus theoreti- 
cally: 


a 


d 



( 8 ) 


There is yet no comparable expression for h . 

w 

The tropospheric contribution to a vertical range 
measurement on the basis of Eqs. (5) and (7) then becomes: 


Ah, = Ah + Ah 
tro tro, tro 

d w 


Ah = 10 
tro , 


-6 _ h d 


0 


.-6 1 


J~N d dh=10 .-N d a d T s \(9) 


Ah 


tro 


w 


o 


N dh = 10“ 6 • i N h 
w 5 w w 
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All the parameters for Eqs. (6) and (9) can be ob- 
tained by least squares procedures from observed zenith 
integrals of N that are obtained from balloon meteorologi- 
cal data. Preliminary values for several geographical lo- 
cations were reported earlier (Refs. 6 and 7). 

In the present study, the approximations made in the 
earlier work have been improved. Revised parameters are 
given for all the former stations and for some new ones. 
Further development of the theory leads to a new method of 
estimating the percent of water vapor in the atmosphere, 
and results are given. Some statistics on zenith range ef- 
fects and some preliminary work on horizontal gradients 
are also presented. 
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2. THEORY 


The theoretical expression for k of Eq. (6), derived 
earlier (Ref. 6), is independent of temperature or its lapse 
rate, and is: 


k 


77. 6 R „„-9 

x 10 

S 


( 10 ) 


where all quantities are in cgs units. The height variation 
of g within the nonionized atmosphere is small and may be 
neglected if a value of g at midatmosphere is used (~ the 
500-millibar level for a sea-level location). Since k is a 
function of g, its theoretical values show a latitude varia- 
tion. 


In order to deduce the value of k from observed 
data, reliable values of dh are needed, to be used in 
Eq. (6). The following improvements have been made in 
the earlier assumptions and computing procedures for de- 
termining Jn^j dh: 

1. Initially, the height integrals of N used the 
"geopotential heights" of the balloon data. These are es- 
sentially determinations of geopotential differences from 
meteorological data. Heights in geopotential meters are 
numerically equivalent to heights in real or geometric 
meters only when g = 980 cm/s^. A height of 10 "geopo- 
tential kilometers" is 10. 036 km at the equator but 9. 983 
km at a pole. The initial "height" integrals were therefore 
an approximation to true or geometric height integrals, and 
although the initial integrals showed a latitude variation 
(Ref. 7), they did not show the theoretical amount. In the 
present work, geometric heights are derived from the 
"geopotential heights" of the balloon data (Ref. 8) and are 
used in the numerical integration. 
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2. The contribution of the unobserved top part of 
the atmosphere can now be approximated better than be- 
fore, using the results of Refs. 6 and 7. 


3. In the earlier numerical integration with 
height, a linear decrease of N d was assumed between ob- 
served points. This would bias the integral slightly too 
high. An approximate correction has been developed on 
the basis of a nominal profile shape. 

Details of these three improvements will be given 
in the Section "Data and Computations. " 

The empirical values of k are based on the real 
atmosphere, which contains a small percent of water 
vapor. If theoretical values are computed from Eq. (10) 
on the basis of the gas constant, R, per gram of dry air, 
a small discrepancy may be expected and may be used to 
provide a measure of the water vapor content of the at- 
mosphere. 

The gas constant, R, per gram of air (observed 
sample) is: 


R 


R 


❖ 


M 


obs 


( 11 ) 


where M Q k g is the mean "molecular weight" of the ob- 
served sample of air and R' 1 ' is the gas constant per mole. 
Combining Eqs. (10) and (11) and solving for M , , we 
obtain: ODS 


M 


obs 


10~ 9 x 77. 6 R* 
gk 


(12) 


Thus a value of M 0 b s is obtained for each set of data that 
provides a value of k. 
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Now let: 

x = mean mole fraction of water vapor in 
the air sample, 

1 - x = mole fraction of dry air, 

M = molecular weight of water = 18.0160 
w 

M = molecular weight of dry air = 28.966 
(Ref. 8). 

Then taking M 0 ^ g as a weighted mean of and M w : 


M x + M,(l - x) = M , 
w d obs 


and solving: 


M - M , 
d obs 

M . - M 
d w 


(13) 


(14) 


Although the zenith integral of can be predicted 
from surface pressure alone with Eq. (6), prediction at 
large zenith angles must also involve the temperature and 
its height variation. The quartic profile, which fits a 
temperature lapse rate of 6. 8°C/km when g = 980 cm/s^, 
is suitable for most regions. For empirical study, the 
height, h^, is written: 



a.T , 
d c 


(15) 


where h, is the "equivalent height" of the model at a sur- 
d 0 

face temperature of 0°C and a^ is its temperature coeffi- 
cient. For the wet component, there is as yet no compara- 
ble expression, and we assume that is a constant for a 
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given data sample. Integrating Eqs. (5) with these values, 
these theoretical integrals are compared with the revised 
integrals obtained from the balloon data in order to get im- 
proved values of h d , a d , and h^ for each data set, by a 


least squares procedure. Prediction of the zenith integrals 
from these values will be reported below; its precision is 
independent of whether or not the shape of the model profile 
matches the observed profile well. This is not true at low 
elevation angles. The prediction errors to be expected at 
low angles have not yet been studied. 


Equations (2), (3), and (4) would make it possible to 
write a set of mathematical expressions to correspond to 
any observed N d profile. The profile would not generally 
be a continuous function but would consist of a sequence of 
curved segments, one segment for each air layer having 
the same lapse rate. The equations would have the form of 
Eq. (2) with a different value of p, for each layer; they 
would be easy to write but not necessarily easy to handle, 
since n would generally not be an integer. The height 
parameter for each layer, however, would be deduced easily 
from the lapse rate in the layer and the temperature at the 
base of the layer. The agreement between Eqs. (3) and (4), 
and the empirical parameters found for Eq. (15), indicate 
that this procedure could give a close approximation to the 
observed N d profile. 
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3. DATA AND COMPUTATIONS 


The data were obtained from meteorological balloon 
ascents and were supplied by the U.S. National Climatic 
Center. The refractivity at each height of observation was 
computed from Eqs. (1), as in Ref. 6. The height integrals 
of N^ and N w were then obtained by numerical integration. 
As before, data from a given balloon were used only if ob- 
servations were complete up to the 30 -millibar level or 
higher for pressure and temperature, and up to the 500- 
millibar level for water vapor. 

The three improvements described below were made 
in computing the numerical integrals. The changes were 
made one at a time so that the individual effect of each could 
be noted. 

1. The geometric heights were computed from the 
"geopotential heights" of the data, using the equations and 
parameters of Ref. 8. The procedure takes into account 
the ellipsoidal shape of the earth, its mean density distri- 
bution, centrifugal acceleration, and the inverse square 
law, but not local gravitational anomalies. The procedure 
is considered satisfactory even up to heights of several 
hundred kilometers. The heights involved in the balloon 
data are all less than 40 km. The resulting geometric, or 
true, height integrals differ from the earlier "geopotential 
height" integrals of Refs. 6 and 7 by latitude -dependent 
amounts up to 8 or 9 mm (out of a total integral of 2. 3 
meters). The change involves an increase near the equator 
and a decrease near the poles. 

2. The new expression for the contribution of the 
unobserved top part of the atmosphere is: 

-3 

Ah . = 2.296 x 10 P, (meters), 

tro(top) top 
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where P^ 0 p is the topmost observed pressure in millibars. 
This is intended to be used at heights not less than 23 km 
(30 -millibar level), and the constant factor is based on a 
nominal value of g in that height region. This change pro- 
duces a small but observable improvement in the precision 
of the least squares fit in getting parameters. 

3. In the earlier numerical integration of with 
height, it was assumed that the decrease of is linear 
between successive observed heights. The resulting inte- 
grals all must be slightly too large. To improve this, it is 
here assumed that decreases as a fourth-degree function 
up to a height of 5 km, and decreases exponentially above 
that (Appendix A). A small curvature correction is computed 
for each layer (height interval) and subtracted from the 
layer integral obtained by the earlier procedure. This im- 
provement decreases the total integral Jn^ dh by approxi- 
mately 2. 0 ± 0. 5 mm. The larger of these curvature cor- 
rections occur, as expected, for profiles having fewer 
observed points. The rms fit of the derived parameters 
is slightly improved. This change has not been applied to 
the integration of N w , where the observed values of N w 
are not always monotonic with height. 

All the data were in 1-year sets (two balloons per 
day, at 0^ and 12^ UT at all stations). Several data sets 
have been added to the list of Ref. 6. The results given 
here cover the complete set of stations used and come from 
the improved computing procedures. 
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4. RESULTS AND DISCUSSION 


COEFFICIENT k FOR EQ. (6), AND INFERRED WATER 
VAPOR CONTENT 

Equation (6) states a linear relation between Jn , dh 
and the surface pressure. The coefficient, k, was deter- 
mined for each 1-year set of data separately. The results 
for all the stations are listed in Table 1. The rms error 
in predicting the height integral from this value of k (also 
listed) is between 1 and 2 mm for each 1-year data set. 
Figure 1 shows representative samples of observed and 
predicted integrals, J N d dh, for comparison. 

Figure 2 shows these observed values of k as a 
function of latitude; it also shows the theoretical k values 
for a dry atmosphere, calculated from Eq. (10), using the 
local value of g computed from geodetic parameters 
(Applied Physics Laboratory geodesy model No. 4.5) for 
the observed local height of the midpressure level. These 
calculated points do not fall on a smooth curve because of 
local field variations and the variety of station heights 
above the geoid (especially the latter). If the midatmo- 
sphere height were everywhere the same (e.g., 5, 6, or 
7 km), and if the earth were an ellipsoid, one of the solid 
curves of Fig. 1 would represent the theoretical value of 
k. In fact, however, the usual height of the half-pressure 
level above a sea-level station is approximately 6 km near 
the equator and 5 km near a pole. If the station is much 
above sea level, the midpressure level is naturally higher 
also, about 1 km higher than would otherwise be expected 
if the station height is 1. 5 km (see Albuquerque and Byrd 
Station). 

The bias between observed and theoretical values 
of k for a dry atmosphere is conspicuous and is markedly 
dependent on latitude. As was discussed above, theory 
shows that this difference may be due to the moisture con- 
tent of the real atmosphere. The average fractional con- 
tent of water vapor molecules in the atmosphere, during 
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• OBSERVED 
O PREDICTED FROM P s 
COLUMBIA, MISSOURI 



DAY OF MONTH (UT) 

Fig. 1 OBSERVED AND PREDICTED VALUES OF /N d dh 
FOR COLUMBIA, MISSOURI, JANUARY AND 
JULY 1967 
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J OBSERVED ± a (o OMITTED FOR VERY CLOSE POINTS) 

A CALCULATED USING g AT OBSERVED ATMOSPHERE CENTROID 
APL GEODESY, DRY AIR 


CALCULATED FOR ELLIPSOIDAL EARTH, Dl FFERENT CENTROIDS 

DRY AIR 



Fig. 2 LATITUDE VARIATION OF TROPOSPHERIC RANGE 
ERROR PARAMETER k (WHERE /N d dh = kP s ) 
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the year, has been computed on this basis for each 1-year 
data set from Eqs. (12) and (14). The results are given in 
Table 2 and plotted in Fig. 3. The values must be consid- 
ered approximate, since any data bias that leads to an 
error in M^g will obviously have a large effect on the de- 
rived value of the water vapor content. The values, how- 
ever, are of reasonable magnitude. Further evidence that 
they are at least qualitatively correct will be shown in a 
later figure. 


PARAMETERS FOR THE TWO-QUARTIC N MODEL 


The revised computation procedures were also used 
to get revised height parameters h^ , a d , and h^. for the 


N model for all stations (Eq. (5)). 
listed in Table 3. 


These parameters are 


The zenith integral, J* N d dh, can be predicted from 
the model and from these parameters, at least as precisely 
as from the surface pressure and the k values, as will be 
seen from the residual errors of prediction in Tables 1 
and 3. 


No detailed study has been made of the latitude ef- 
fects and water vapor effects that might be expected in 
these "dry" parameters. However, some nominal theo- 
retical values can be given for comparison with the mean 
values of Table 3. 

Starting with Eq. (4), the use of a quartic N d model 
implies a temperature lapse rate of 6. 814°C/km if g = 978 
cm /s^ (midlatitude, midatmosphere value). From this, 

Eq. (3) gives h H = 273. 16/6. 814, or 40. 088 km, which is 
a 0 

in good agreement with the values of Table 3. From Eq. 
(8), a d = 0. 14675 km/°C, whereas the mean of Table 3 is 
0. 14872. Theoretical values based on dry air and nominal 
values of g can thus provide a good approximation to the 
observed parameters. 
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Table 2 


Water Vapor Content, x, from k Values 
(Fraction of Total Air Molecules) 


Station 

Year 

Height of 

Atmosphere Centroid 
(midpressure) 
(km) 

g 

at Centroid 
(cm/s^) 

Water Vapor 
Fraction, x 

Weather Ship E 

1963 

5. 83 

977. 95 

0.0045 

Weather Ship E 

1965 

5. 83 

977. 95 

0.0042 

Weather Ship E 

1967 

5. 83 

977. 95 

0.0040 

Ascension Island 

1967 

5. 85 

976. 36 

0.0106 

Caribou, Maine 

1967 

5.52 

979. 06 

0.0031 

Washington, D. C. 

1967 

5. 68 

978. 30 

0.0040 

St. Cloud, Minn. 

1967 

5.57 

978. 97 

0.0035 

Columbia, Mo. 

1967 

5. 70 

978. 33 

0.0044 

Albuquerque, 
New Mexico 

1967 

7.05 

977. 56 

0.0026 

El Paso, Texas 

1967 

6.90 

977. 37 

0.0041 

Vandenberg 
AFB, Cal. 

1967 

5. 82 

977. 90 

0.0036 

Pago Pago, 
Samoa 

1967 

5. 86 

976. 47 

0.0076 

Wake Island 

1963 

5. 80 

976. 87 

0.0068 

Wake Island 

1965 

5. 80 

976. 87 

0.0070 

Wake Island 

1967 

5. 80 

976. 87 

0.0071 

Majuro Island 

1967 

5. 87 

976. 28 

0.0090 

Point Barrow, 
Alaska 

1967 

5. 32 

981.04 

0.0034 

Byrd Station, 
Antarctica 

1967 

6.63 

981. 10 

O'. 0020 
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LATITUDE (degrees) 

Fig. 3 MEAN WATER VAPOR CONTENT OF ATMOSPHERE, 
VERSUS LATITUDE 
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Prediction of the wet component, f N w dh, is much 
poorer and leaves a residual error usually of a few centi- 
meters instead of a millimeter or two. This problem will 
be discussed later. 


STATISTICS ON ZENITH INTEGRALS 

Statistical values of the zenith integral are needed 
for range correction when no surface data are available 
(e.g.j radar altimeter over ocean areas), and they are 
also of climatological interest. Table 4 gives 1-year mean 
values and standard deviations for the dry and wet com- 
ponent integrals and the totals for all the stations studied, 
based on the revised computations described above. 

The rms error of using the annual mean to repre- 
sent the height error at a specific time varies for different 
stations from a few millimeters to 2 cm for the dry com- 
ponent, and from 3 to 6 cm for the wet component, except 
at the Antarctica station (where it is less than 1 cm). The 
larger percent variations from the mean, in both compo- 
nents, are generally found for the higher -latitude stations, 
where seasonal changes are more marked. Statistical 
prediction can probably be improved by including a sea- 
sonal effect, but such parameters have not yet been ob- 
tained. 


If there is no correlation between dry and wet com- 
ponent integrals on a day-by-day basis, the variance of 
the total integral should be equal to the sum of the vari- 
ances of the components. On this basis, no conspicuous 
evidence of correlation between dry and wet integrals is 
seen in the present statistics; here, the variance of the 
total at one location is not markedly different from the 
sum of the component variances, and the deviation is •not 
always in the same direction. However, this approach may 
not be suitable for detecting small effects. Some relations 
between the dry and wet parts have in fact been observed 
but have not yet been fully treated theoretically. 
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Table 4 


Annual Statistics on Jn dh and Components 




J* dh (meters) 

J*N dh (meters) 

Total Jn dh (meters) 



Mean 

<7 

Mean 

a 

Mean 

CT 

Weather Ship E 

1963 

2. 32546 

0. 01521 

0.21370 

0. 04996 

2.53916 

0.05538 

Weather Ship E 

1965 

2. 32202 

0.01694 

0. 18578 

0.05174 

2. 50779 

0.05888 

Weather Ship E 

1967 

2. 32843 

0.01493 

0. 16647 

0. 05038 

2. 49490 

0.05151 

Ascension Island 

1967 

2. 30176 

0. 00563 

0. 17516 

0.03105 

2. 47692 

0.02889 

Caribou, Maine 

1967 

2. 26051 

0. 01914 

0.08675 

0. 05839 

2. 34727 

0.05918 

Washington, D.C. 

1967 

2, 29975 

0. 01545 

0. 11671 

0. 06485 

2. 41646 

0.06366 

St. Cloud, Minn. 

1967 

2. 22660 

0. 01800 

0.08252 

0. 05224 

2. 30912 

0.05140 

Columbia, Mo. 

1967 

2. 25361 

0.01440 

0. 10922 

0.06167 

2. 36283 

0.05978 

Albuquerque, 
New Mexico 

1967 

1.91237 

0. 01132 

0.06752 

0.04321 

i 

1. 979 89 

0.04797 

El Paso, Texas 

1967 

2.01287 

0.00986 

0.08923 

0.05047 

2. 10211 

0.05252 

Vandenberg 
AFB, Cal. 

1967 

2.29101 

0.00882 

0.09298 

0.03329 

2.38399 

0.03202 

Pago Pago, 
Samoa 

1967 

2. 31367 

0.00492 

0.25088 

0.05089 

2. 56455 

0.05101 

Wake Island 

1963 

2. 31700 

0. 00454 

0.23929 

0.05230 

2. 55628 

0.05107 

Wake Island 

1965 

2. 31980 

0.00431 

0. 19655 

0. 04698 

2. 51635 

0.04651 

Wake Island 

1967 

2. 31756 

0.00534 

0.21938 

0. 05774 

2.53694 

0.05631 

Majuro Island 

1967 

2. 31261 

0.00354 

0. 27775 

0. 05354 

2. 59036 

0.05417 

Point Barrow, 
Alaska 

1967 

2. 30310 

0. 02024 

0.04712 

0.02862 

2. 35022 

0.03037 

Byrd Station, 
Antarctica 

1967 

1. 85214 

0.02115 

0.01520 

0.00696 

1. 86734 

0.02461 
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OBSERVED RELATIONS BETWEEN DRY AND WET 
COMPONENTS 

The mean percent of water vapor molecules in the 
atmosphere at each station was deduced above from the 
dry component integral, jN d dh (Table 2). The mean wet 
component integrals, Jn w dh, for all stations are listed in 
Table 4. These two quantities are independently measured 
and are not measurements of the same thing; still, they 
both involve the water vapor content, and it is of interest 
to compare them. Figure 4 shows the comparison. The 
wet component integrals show a clear, though not precise, 
relation to the water vapor content as determined from the 
dry integrals. This corroborates the assumption made 
above that the systematic discrepancy between theoretical 
and observed values of the parameter k (Eq. (10)) is 
caused by the presence of water vapor in the atmosphere. 

Another bit of corroborating evidence is found in 
the residual errors when jN d dh is predicted from surface 
pressure using Eq. (6) and the k values of Table 1. The 
observed and predicted values of Jn^ dh were shown in 
Fig. 1 for Columbia, Missouri, for one summer and one 
winter month. Each observed point pertains to one balloon 
flight. The prediction errors are small but systematic; 
predicted values are too high in January and too low in 
July. These residual errors on a much-enlarged scale 
are plotted in Fig. 5 as a function of observed water vapor 
pressure at the surface at the time of the balloon flight. 

The data scatter appears large on this large scale, but the 
trend is obvious. Data from all stations show the effect, 
but it is most observable where the seasonal variation of 
water vapor pressure is large. 

The prediction of Jn^ dh can be improved on the 
basis of these findings. Empirically, we may rewrite 
Eq. (6) in the form; 



i.5 


WATER VAPOR CONTENT, DEDUCED FROM /N d dh (percent) 

Fig. 4 COMPARISON OF/N w dh THROUGH ATMOSPHERE, 
WITH WATER VAPOR CONTENT DEDUCED 
FROM /N^dh 
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Fig. 5 WATER VAPOR DEPENDENCE OF RESIDUAL ERRORS 
OF PREDICTED/N d dh 
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where kj and k 2 are constants, and e gur is the observed 
value of e at the surface. Preliminary values of kj and 
k 2 have been obtained from the balloon data, again by a 
least squares procedure. Equation (16) improves the pre- 
diction of jN d dh. Its use removes the seasonal bias from 
the prediction and results in smaller residual errors: gen- 
erally an rms value of 1. 5 mm or less at the zenith. 

Values of kj are slightly smaller than those of k for the 
same location, and k 2 has approximately one -tenth the 
value of kj. 

Equation (16) has the disadvantage that the measure- 
ment of e gur is much less precise than that of P g , and this 
lack of precision degrades the computed value of kj as 
well as k 2 . Further work is needed here, and the prelimi- 
nary values of kj and k 2 will therefore not be listed. The 
theoretical aspect of Eq. (16) should be related to the ap- 
proximations that were originally made in defining N d 
(Ref. 1). This question has not been studied thoroughly 
and will not be discussed further at this time. 


HORIZONTAL GRADIENTS IN THE ATMOSPHERE 

Up to this point, only vertical gradients of meteoro- 
logical quantities have been considered; it has been as- 
sumed that horizontal gradients are negligible. Several 
figures will be presented here to give a qualitative idea of 
the relation between simultaneously observed zenith inte- 
grals at locations a few hundred miles apart. 

Figure 6 shows the observed and calculated values 
of the dry zenith integral as a function of time for the 
months of January and July 1967 at Albuquerque, New 
Mexico, and El Paso, Texas, 200 miles to the south. For 
each month, the data for the two stations are plotted on 
the same scale, with no zero shift. The bias between 
Albuquerque and El Paso data is due to the difference in 
altitude of the stations. Albuquerque is 400 meters the 
higher of the two, with correspondingly less air above it: 
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DAY OF MONTH 

Fig. 6 COMPARISON OF VALUES OF /N d dh AT EL PASO 
AND ALBUQUERQUE, WINTER AND SUMMER 
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the zenith range effect (dry component) at Albuquerque is 
10 cm smaller than that at El Paso. The similarity of the 
pattern for the two stations is very marked at both seasons, 
the indication being that the horizontal variations are not 
very large over a few tens of miles. The consistency of 
the pattern for the two stations speaks well for the preci- 
sion of the data. Each observed point was obtained from 
a different balloon flight, each with its own measuring 
instruments. 

A diurnal variation is seen, especially in the sum- 
mer samples (both observed and predicted). It is larger 
than the expected atmospheric tidal effect and is not yet 
explained. 

Figure 7 shows corresponding observed data for 
January 1967 for three other U. S. stations: St. Cloud, 
Minnesota; Washington, D. C. (Dulles Airport); and 
Caribou, Maine. With Columbia (Fig. 1), these form the 
vertices of a quadrilateral that is roughly trapezoidal, and 
their separations are from 500 to about 1200 miles. Co- 
lumbia is some 800 miles from Albuquerque. Over these 
distances, there is not the same close correlation as in 
Fi-S* 0, but marked similarities can often be seen even at 
distances greater than 500 miles, sometimes with a time 
delay (not always in the same direction). The deep mini- 
mum that appeared at Washington on 2 8 January and at 
Caribou a day later is a conspicuous example. 

Figure 8 shows the same three stations as Fig. 7 
but for the month of July, and Figs. 1 and 6 again supply 
further data. Correlations between stations are very 
noticeable in Fig. 8. The zenith integral (dry) was more 
variable in the January sample than in the July sample at 
all these stations. On a large-scale basis, it appears 
that the horizontal gradients are greater in winter, but 
this conclusion cannot be fully justified without more work 
on the apparent diurnal variation. 


Less consistency is found between zenith integrals 
of the wet component at different stations, but some pattern 
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2.30 


2.25 


2.20 



JULY 1967 (days) 

Fig. 8 COMPARISON OF SUMMER SAMPLES OF /N d dh 
AT THREE LOCATIONS IN THE U.S. 
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similarities can be observed. Figures 9 and 10 show a 
winter sample and a summer sample, respectively, for 
the St. Cloud, Caribou, and Washington stations. The 
difficulty of measuring water vapor content is a serious 
problem here, and improved methods are needed. It is 
impossible to be sure how much of the variation in these 
data is due to atmospheric variation and how much to the 
measurement. It would be difficult to get any reliable in- 
formation about horizontal gradients from these data. 

As expected, the wet zenith integral is considerably 
greater in summer than in winter at either station, and is 
greater on the average in the Washington area than in 
Maine at either season. There must be a systematic (but 
small) latitude gradient of the mean value of J*N W dh,as 
well as systematic climatic gradients owing to geography 
(e.g. , land versus ocean areas). 

In recent work at the National Severe Storms La- 
boratory (Refs. 9 and 10), horizontal gradients of pres- 
sure, temperature, and water vapor pressure have been 
examined by the simultaneous launch of meteorological 
balloons about 15 nmi apart. The statistical results indi- 
cate that the pressure difference between stations at this 
separation is too small to be measured with present instru- 
ments. Detection of a temperature difference is marginal, 
but the statistics indicate that differences larger than 1° or 
2°C are not commonly observed. Water vapor statistics 
are more variable and also more uncertain because of 
low instrumentation accuracy (Ref. 10). 

Clearly it will be impossible to estimate Jn w dh 
with accuracy comparable to the J N d dh prediction until 
the wet integral can itself be better measured than at pres- 
ent. It is not known how much of the variability now ob- 
served in the values of Jn w dh can reasonably be attributed 
to the present measuring errors. This question will be 
studied further. 

Additional instrumentation could help to solve this 
problem. If a refractometer were added to measure the 
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total refractivity, N, along with the usual meteorological 
quantities, then the wet component, N w , could be found 
from the difference, N - N d , and compared with the value 
from humidity measurement. Experimentation with this 
redundant instrumentation at a surface site could deter- 
mine whether its use would be justified in balloon flights. 
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APPENDIX A 

NUMERICAL INTEGRATION ALONG A CURVE: 
CORRECTION TO TRAPEZOIDAL RULE 



Numerical integration of N with height from a set 
of observed points is most simply accomplished by assum- 
ing a linear change of N with height between successive 
points. For a function whose curvature is practically al- 
ways in the same direction (e. g. , N d versus height), the 
integral so computed will be slightly biased. The correc- 
tion described here is not used in integrating N , which has 
a less consistent profile. 

The separations. Ah, between observed points in 
these profiles varied widely (from a few meters to 2 km), 
and the theoretically correct functional form of the N d 
curves varied also. Expressions for a correction were de- 
rived as given below, based (1) on a fourth-degree N d pro- 
file (temperature lapse rate 6. 8°'C/km), and (2) on an ex- 
ponential profile (zero lapse rate). The first correction was 
applied to the lower and the second to the upper part of the 
atmosphere. The area under a chord joining two observed 
points is: 

,h 2 N + N 

chord: J N dh = - Ah . (A-l) 
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For case (1), let both points lie on a fourth-degree function 
given by (Ref. 5): 


N 


N. 


< h d - h i> 



(A-2) 


the vertex being at h d . Let l = h - hj. Then using values 
from Eq. (A-2) in Eq. (A-l), it can be shown that: 


chord: 


/ 


N 


N dh = 


* 1 
- 2 i 


V( 


4 q o 9 1 

l Ah - 21 Ah + 31 Ah 


4 — 5\ 

- + *£-)■ 


(A -3) 


Integrating the fourth-degree curve of Eq. (A-2) be- 
tween the same limits yields: 


r 2 N i/ 

/ Ndh=-(i 


quartic: / N dh 

h l 


Ah - 21 3 Ah + 2£ 2 Ah 


- i Ah 4 + 


)• 


(A -4) 


The difference (chord - quartic) is the correction, and is: 



Since Ah is always small in comparison with i, the first 
term is by far the most significant. 


For case (2), let two points lie on an exponential 
curve given by: 


N = Nl e- c(h - hl) . 


(A-6) 
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Then 

N 2 = N l e " CAh * (A-7) 

and using Eq. (A-l), the integral along the chord is: 
r 2 N 1 Ah 

chord: J Ndh = — (1 + e‘ C ) Ah. (A-8) 


we get: 


Integrating Eq. (A-6) between the same two points. 


2 N 

jxponential: J' Ndh = (1 - e C 


(A -9) 


The difference between Eqs. (A-8) and (A -9) then gives the 
correction. Expanding e -c in the usual infinite series 
in both equations and subtracting, the difference (linear - 
exponential) becomes: 

**2 N 

A j N dh = ~y c 2 Ah 3 || - | c Ah+ ...j .(A-10) 
h ! 

In order to evaluate and compare the quartic and the 
exponential corrections for curvature for a given interval. 
Ah, values for the parameters i in Eq. (A -5) and c in 
Eq (A-10) are needed. It follows from the results of Ref. 6 
that in order to give the correct integral, Jn^ dh, through 
the atmosphere above any height, hj, the parameter i for 
a quartic profile is: 


2 




(A- 11) 
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where T. is the Celsius temperature at hj, j0q = 40 km 
(the equivalent height of the quartic model at 0°C), and a = 
0. 147 km/°C. Thus i can be computed for any data point. 


For the quartic model of Eq. (A -2): 


h, + i 


/ 


N dh 


Nji 
~ ' 


(A- 12) 


For the exponential model of Eq. (A-6): 



N dh 



(A - 13) 


If the exponential model is to yield the same integral above 
hj as the quartic model, then: 


c 



(A- 14) 


The corrections of Eqs. (A -5) and (A- 10) can thus be found. 
Since they are small and are approximations at best, it is 
not necessary to compute them precisely. It was estimated 
from observed profile data that the corrections for any in- 
terval, Ah, can be approximated adequately by the values: 


(linear -quartic curve) A 



N dh 



(A- 15) 


(linear -exponential curve) 



N dh = 1. 9 


-^Ah 3 . (A- 16) 

t 


Either correction is to be subtracted from the corre- 
sponding integral obtained from the trapezoidal rule. 


- 38 - 



THE JOHNS HOPKINS UNIVERSITY 

APPLIED PHYSICS LABORATORY 

Silver Spring. Maryland 


REFERENCES 


1* E. K. Smith, Jr. and S. Weintraub, "The Constants 
in the Equation for Atmospheric Refractive Index at 
Radio Frequencies, " Proc. IRE , Vol. 41, No. 8, 
August 1953, p. 1035. 

2. B. Edl4n, "The Dispersion of Standard Air, " J. Opt. 
Soc. Am. , Vol. 43, No. 5, May 1953, p. 339. 

3. B. R. Bean and E. J. Dutton, "CRPL Exponential 
Reference Atmosphere, " U. S. National Bureau of 
Standards Monograph No. 4, 1959. 

4. B. R. Bean, "Concerning the Bi-exponential Nature 
of the Tropospheric Radio Refractive Index, " Beitr . 
Fhys. Atmos. , Vol. 34, 1961, p. 81. 

5. H. S. Hopfield, "Two-Quartic Tropospheric Refrac- 
tivity Profile for Correcting Satellite Data, " J. 
Geophys. Res ., Vol. 74, No. 18, 20 August 1969, 

p. 4487. 

6. H. S. Hopfield, "Tropospheric Effect on Electro- 
magnetically Measured Range: Prediction from 
Surface Weather Data, " Radio Sci . , Vol. 6, No. 3, 
March 1971, p. 357. 

7. H. S. Hopfield, Tropospheric Range Effect at the 
Zenith, JHU/APLTG 1168, June 1971. 

8. Smithsonian Meteorological Tables , 6th Rev. Ed. , 

R. J. List, Ed., 1963. 


- 39 - 



THE JOHNS HOPKINS UNIVERSITY 

APPLIED PHYSICS LABORATORY 
Silver Spring. Maryland 


REFERENCES (Cont'd) 


9. S. L. Barnes, J. H. Henderson, and R. J. Ketchum, 
Rawinsonde Observation and Processing Techniques 
at the National Severe Storms Laboratory , NOAA 
ERL NSSL-53, April 1971. 


10 . 


S. L. Barnes (private communication), 1971. 



